Energy and momentum deposited into a QCD medium by a jet shower 
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Hard partons moving through a dense QCD medium lose energy by radiative emissions and 
elastic scatterings. Deposition of the radiative contribution into the medium requires rescattering 
of the radiated gluons. We compute the total energy loss and its deposition into the medium 
self-consistently within the same formalism, assuming perturbative interaction between probe and 
medium. The same transport coefficients that control energy loss of the hard parton determine 
how the energy is deposited into the medium; this allows a parameter free calculation of the latter 
once the former have been computed or extracted from experimental energy loss data. We compute 
them for a perturbative medium in hard thermal loop (HTL) approximation. Assuming that the 
deposited energy-momentum is equilibrated after a short relaxation time, we compute the medium's 
hydrodynamical response and obtain a conical pattern that is strongly enhanced by showering. 

PACS numbers: 12.38.Mh, ll.lO.Wx, 25.75.Dw 



Jet quenching (the modification of hard jets in dense 
media) is one of the most studied discoveries at the Rela- 
tivistic Heavy-Ion Colhder (RHIC) fl. It is expected to 
play a key role in the study of the quark-gluon plasma 
(QGP) produced in heavy-ion collisions at the Large 
Hadron Collider (LHC). Numerous experiments Q have 
established the suppression of hadrons with high trans- 
verse momenta; others indicate that the lost energy man- 
ifests itself as conical flow in the soft sector 0] . 

Calculations of jet modification tend to focus on one 
of two separate questions: the modification of the final 
hadron distribution from the hard parton due to its en- 
ergy loss, or the response of the medium to the energy 
deposited. Numerous studies of the former, based on per- 
turbative QCD (pQCD), have yielded near-rigorous mea- 
sures of the two non-perturbative transport coefficients 
q = ^IjY^ and e = ^ which codify the transverse (to the jet 
axis) momentum diffusion and longitudinal drag experi- 
enced by a fast parton Q . Computations of the medium 
response consist of two parts: an ansatz for the space- 
time profile of the energy-momentum deposition, and a 
calculation of the dynamical response to this "source" of 
excess energy and momentum. Based on its success at 
RHIC, ideal fluid dynamics has been used to compute 
this medium response assuming that the energy lost 
by the jet is entirely deposited into the medium at a con- 
stant rate and thermalizes instantaneously. 

So far there exists no first principles calculation of 
the magnitude and space-time profile of the energy- 
momentum deposited in a medium by a hard parton that 
can be considered on par with the pQCD energy loss 
calculations A noteworthy attempt to calculate the 
deposition profile in pQCD is the semi-phenomenological 
approach of Neufeld and Miiller who use the differen- 
tial single gluon emission spectrum of Ref. 0] and inter- 
pret this as the rate of gluon emission in the medium. A 
non-diffusive Fokker-Planck equation is then motivated 
to compute how this distribution changes due to elastic 



energy loss of the emitted gluons. As anticipated in [9|, 
they find that not all of the energy lost to gluon radia- 
tion is deposited in the medium. However, since the un- 
derlying formalism lacks information about virtuality 
evolution, this calculation docs not include gluon multi- 
plication by showering, i.e. the splitting of a radiated 
gluon into two lower virtuality gluons. The transverse 
momentum deposition thus cannot be computed and, due 
to the strict eikonal limit used in Q, the parent parton 
docs not lose energy after radiation. We present a new 
formalism in which the radiative and elastic energy loss 
of the fast parton, its virtuality evolution by radiation, 
the showering and multiplication of the radiated gluons, 
and the energy deposited by them in the medium are all 
calculated consistently in the same approach. 

Hard jets in vacuum or in heavy ion collisions are pro- 
duced with considerable virtuality. As the jets proceed 
through vacuum or medium, this virtuality is lost by se- 
quential radiative emissions. The effect of this perturba- 
tive shower on the non-perturbative hadronization pro- 
cess is computed using DGLAP evolution equations 
for the fragmentation function. These equations express 
the radiation of multiple partons, which hadronize inde- 
pendently, via an evolution in virtuality of the parent 
parton. In a medium, one can derive analogous equa- 
tions where the gluon radiation probability is modified 
by the scattering of hard parton and emitted gluons off 
medium constituents. These are referred to as "medium 
modified evolution equations" . In addition to stimulating 
gluon emission, the scattering of the hard parton causes 
it to lose forward light-cone momentum by elastic ex- 



changes with the medium [llj, [12]. At the same time 



the parton gains transverse momentum from the medium 
[Tst and imparts to it an equal amount in return. In 
an arbitrary medium these effects are encoded in two 
non-perturbative transport coefficients, q and e, defined 
in terms of in-medium gluon field correlation functions 
13, The medium modification of the standard vac- 
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uum evolution depends on these transport coefficients. 

In-medium evolution equations where the medium 
modified fragmentation function (MMFF) is affected only 
by q were derived in 1J| . We point out that the same pro- 
cesses can be used to compute the amplification of the en- 
ergy deposited through multiple radiations stimulated by 
transverse broadening. Formally, this can be computed 
by replacing the operator expression for the fragmenta- 
tion function with that for the energy deposited; this is 
identical to e. Using this calculation of the modified e 
as the energy deposited through all elastic scatterings of 
the shower places it on the same footing as energy loss. 
The diagrams involved and the resulting expressions for 
the in-medium splitting functions (IMSF) are identical. 
The solution of the in-medium evolution equation for e 
no longer represents the elastic energy loss by one parton, 
but rather the energy deposited by the jet shower. 

Imagine a hard quark or gluon with large light cone 
momentum q~ (and thus energy E = q~ /^/2) and virtu- 
ality < fjL entering a medium of fixed length L held at a 
constant temperature T. Let us assume that the rate of 
energy deposition by this jet in the medium as a function 
of length C, denoted as ^^^^{EX, IJ^), is known (i.e., can 
be calculated or measured). Note that both the deposited 
energy AE and C, are actually the light-cone quantities 
Ag" and For brevity we refer to these simply as de- 
posited energy and distance travelled. Given the above 
function, the total energy deposited by a jet originating 
at location C,i and propagating to C,f is given as 



AE{E,^?f^^. 



dAE 

C,. dC 



(i^,C,M^)''-™ 



(C/~0)e, (1) 



where the last approximate equality is solely for the case 
of a single parton propagating without radiation. 

If the scale ^ is much larger than Aqcd, the change 
with virtuality in the partonic shower pattern may be 
calculated perturbatively: a leading quark at the higher 
virtuality may split into a quark and a gluon with lower 
virtuality, and similarly for a gluon. As a result, there 
is change in the energy deposited in the medium due to 
the increase of the number of partons depositing energy. 
Using the IMSF from IJ], the change in the energy de- 



position by a quark with energy E from Q to Cf due to 
the increase in virtuality fi can be expressed as 



f Q U 



dAE„{E,fi^f/ aJn^)f^ f^f 

;in(/) " - ^/^.//cp,^..(.>c,.^i.)(2) 



AE,{E,^l%+AEg{yE,^l^)lf+AEg{il-~y)E,^l 
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Here the first term in square brackets represents the en- 
ergy deposited by a quark with energy E and virtuality 
/i^, from the initial location Q to the intermediate loca- 
tion C; the second and third terms represent the energy 
deposited by the quark and the emitted gluon with re- 
duced energies yE and (1 — y)E, respectively, from the 



intermediate location C to the final location (f. In Eq. ([2]) 
the quark IMSF Pq^gg is given as [H, [H] 



P 



q^qg 



qCp 1+y^ 
27r/i2 l—y 



2 cos 



K'^Eyil-y) 



(3) 



The increase in the energy deposited due to the split- 
ting of the parton is reduced by the virtual correction 
which restores unitarity to the evolution equations. The 
effect of such corrections on Eq. ([2]) is incorporated by 
subtracting from it the virtual term 



V= 



2tt 



AEg{E,fi^)l^J^dyJ%CPq- 



,(2/,C,M',i?)-(4) 



Along with the energy deposition from a quark jet one 
has to evolve the one from a gluon jet of virtuality < fi, 
using a similar evolution equation that includes the split- 
ting of a gluon into two gluons or a qq pair. Similar to 
the MMFFs, one solves a coupled set of evolution equa- 
tions for AEq{E,ijF')'}.[ and A£'g(£', ^^^C/ ^^^-^ of which 
arc functions of three variables i?, C^, at the scale fj? . 

The evolution equations (|2|3|4p for a quark jet and the 
coupled equations for gluon jets are motivated by exist- 
ing rigorous derivations of the medium modification of 
the fragmentation functions due to gluon emission [l^, 
the accumulation of transverse momentum [isj and longi- 
tudinal drag by propagating hard partons in a QCD 
medium, and the effect of such accumulated momentum 
on radiative processes [l3|. The IMSF ^ accounts for 
interference between diagrams where the gluon is emit- 
ted at the origin or at the location C,. In propagating up 
to C the quark loses a fraction of its energy; while this is 
included in the total energy deposited, its effect on the 
interference pattern in Eq. ([3]) is ignored; this is justified 
in the eikonal limit for the propagating parton. Yet an- 
other approximation is the neglect of the energy lost by 
the radiated (reabsorbed) gluon in the virtual correction. 
Since the radiated gluon in the virtual correction exists 
in only one amplitude, with a single parton in the com- 
plex conjugate, its energy loss is balanced by the quark 
propagating in the loop. 

In the eikonal approximation, the hard jet loses light 
cone momentum and remains close to on-shell, thus the z- 
component of the deposited light-cone momentum is ap- 
proximately equal to the energy deposited {Ap^ — AE). 
Note that the negative light cone momentum [Aq^) is 
not conjugate to C~ and thus it is not inconsistent to 
compute the C,~ dependence of the Ag^ deposited. The 
remaining two components that may be computed are 
the transverse momentum deposited by the jet as a func- 
tion of Q~ . This can again be directly estimated from a 
pQCD calculation: A parton traversing a medium gains 
transverse momentum squared with length as 



<f d{pl)iE,t,^) 

Ci dC 



1 parton 



(C/-C:)«-(5) 
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By momentum conservation this equals the p\ deposited 
in the medium by the same parton. 

For a hard virtual quark the total transverse momen- 
tum deposited increases due to parton splitting. This 
can be calculated using an equation similar to that for 
light-cone momentum deposition. For a quark with en- 
ergy E and virtuality traversing a medium from Qi 
to (J/, the change of the transverse momentum deposited 
with virtuality is obtained as 

rfln(^^) ^ ^^ldyj^ dCP,^,,{yX.l^ ,E) (6) 

The splitting function here is identical to that in Eq. 
and the meaning of the three terms in the bracket is 
analogous to Eq. Further, one must include a virtual 
correction and couple Eq. ([SI) to a similar equation for the 
p\ deposited by a virtual gluon. 

Using Eqs. (|2l6p (along with the coupled ones for gluon 
jets), we can compute the 3-momentum Ag^, de- 
posited by a hard virtual parton, disintegrating into a 
shower of partons, in a dense medium as a function of 
the length C,^ traversed. Similar to the case of in-medium 
evolution equations for the MMFF these equations 
require an initial condition. For the case of the MMFF, 
the only possible choice was to insist that the part of 
the jet with virtuality below a minimum exited the 
medium and use the known vacuum FF at that scale as an 
input. For the deposited part of the energy-momentum 
we here assume that the medium is weakly coupled, thus 
when the virtuality of the parton is fio ~ 4T the de- 
posited energy and p\ can be obtained from the expres- 
sions for e and q in an HTL plasma [ll[: 
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dAE{no,E)/dC 



CB.as{fJ.l)m% log 



{AET/mjj 



2,^1/4 



d{pl){iio,E)/d(: = CRas[^ll)Tml\og[AET/ml]. (7) 

Here m/j is the Debye screening length and Cn is the 
representation specific Casimir. The integrated energy 
and p\ deposited from Eq. ([7]), as a function of the length 
traveled, is plotted for gluons (circled) and quarks as solid 
lines in Figs. [1] and [51 For a consistent description, we 
impose that partons with an energy E < 4T become part 
of the thermal medium. This condition is maintained 
through the evolution equations. 

Using Eq. ([7|) as input, we may calculate the increase in 
the energy and p\ deposition in the medium as a func- 
tion of C, for initially highly virtual hard partons that 
evolve into a radiative shower. Starting from the scale 
of /io = AT (in all calculations we pick T = 300 MeV and 
a partonic plasma with 3 quark flavors) we evolve up to 
an initial scale /i = E/2. These are plotted as dashed 
lines in Figs. [T] and [21 One notes immediately that both 
quantities increase as we evolve up in virtuality. For com- 
parison, we also estimate the total energy lost by the 
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coUisional deposit 

coU.+rad. deposit 

■ - ■ loss, coU.+rad.+flavor 




FIG. 1: (Color online) Dash-dotted: Total energy lost by a 
hard gluon (circled) or quark by radiative, elastic and flavor 
changing processes. Solid: Energy deposited in the medium 
by a hard parton which does not radiate. Dashed: The same 
for a virtual parton devolving into a partonic shower. 
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collisional deposit 

coll.+rad. deposit 




FIG. 2: (Color online) The {p\) deposited by a hard gluon 
(circled) or quark without radiative emission (solid) and with 
a full radiative shower (dashed). 



hard parton due to elastic, radiative inelastic and flavor 
changing interactions (dash-dotted lines in Fig. [J) . The 
last type of energy loss refers to the case where a quark 
splits with the gluon carrying a larger fraction of the 
momentum, or a gluon splits into a quark-antiquark; in 
this case we assume that the entire energy of that parent 
parton has been lost. This leads to a somewhat artificial 
enhancement of the total energy loss. 
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As an illustration of the effect of this energy- 
momentum deposition in the medium, we compute its 
hydrodynamic response to the following source term: 



dC 



dC 



5\r^)5{t-z). (8) 



In this first attempt we ignore the transverse momentum 
contribution to the source current. Following Refs. [isj . 
we assume that the energy deposited is a small pertur- 
bation and solve for the linear response of the medium: 



ST'"'- 



0, 



(9) 



Tq'^ is the unperturbed energy-momentum tensor of a 
homogeneous and static partonic medium in equilibrium. 
The small excess 6T'^'^ is decomposed as 



5T°" = 6e, ST"'=g\ 



(10) 



de is the excess energy density, g is the momentum cur- 
rent density and = is the sound attenuation length. 
For the specific shear viscosity we took f = ^ • We de- 
lay the response to the source J by a time T^-d ~ to 
account for thermalization of the deposited energy. 




FIG. 3: (Color online) The linear fluid dynamical response to 
the energy deposited by a single parton (left) or by a parton- 
initiated shower (right), when the parton is a quark (top) or 
a gluon (bottom). Note the different vertical scales. 

In Fig. [3] we show the azimuthal projection of the en- 
ergy density \x\6e ai t = 5fm/c after the parton is created, 
for a single non-radiating parton (left) and a parton- 
initiated jet shower (right). A gluon (bottom row) de- 
posits more energy than a quark (top row), due to its 
larger color factor that enters both in the clastic energy 
loss and shower production rate. One immediately notes 
that, while the basic Mach cone structure is not changed, 
showering leads to an enhancement by a factor of 3 in the 



overall magnitude of the response. For quark jets our re- 
sults are qualitatively similar to Ref. 0]. 

In this Letter, we have presented a consistent pQCD 
based calculation of the light-cone and transverse mo- 
mentum {Aq~,p^) deposited by a jet in a medium, as 
a function of distance traversed. Assuming a short ther- 
malization time for the deposited energy we also com- 
puted the hydrodynamic response. The pQCD shower 
has the effect of a large part of the energy being de- 
posited later in the history of the jet (3] which tends to 
enhance the Mach cone like structure formed. 

We thank B. Miillcr and R. B. Ncufcld for helpful dis- 
cussions. This work was supported by the U.S. Depart- 
ment of Energy under grant DE-FG02-01ER41190. 



[1] 

[2] 

[3] 

[4] 
[5] 

[6] 



[7] 

[8] 

[9] 
[10] 

[11] 

[12] 
[13] 

[14] 
[15] 
[16] 

[17] 

[18] 



I. Arsene et al, Nucl. Phys. A 757, 1 (2005); B. B. Back 
et ai, ibid. 757, 28 (2005); J. Adams et ai, ibid. 757, 
102 (2005); K. Adcox et al, ibid. 757, 184 (2005). 
K. Adcox et ai, Phys. Rev. Lett. 88, 022301 (2002); 
C. Adler et al., Phys. Rev. Lett. 89, 202301 (2002). 

A. Adare et ai, Phys. Rev. Lett. 98, 232302 (2007); 

B. I. Abelev et ai, Phys. Rev. Lett. 102, 052302 (2009). 
A. Majumder, J. Phys. G 34, S377 (2007). 

A. K. Chaudhuri and U. Heinz, Phys. Rev. Lett. 97, 
062301 (2006); T. Renk and J. Ruppert, Phys . Lett. B 
646, 19 (2007); B. Betz, et al. larXiv:0812.4401l [nucl-th]. 
For an alternate approach that assumes that the hard 
parton couples strongly to the medium see P. M. Chesler 
and L. G. Yaffe, Phys. Rev. Lett. 99, 152001 (2007); 
S. S. Gubser, S. S. Pufu and A. Yarom, ibid. 100, 012301 
(2008). 

R. B. Neufeld and B. MuUer, Phys. Rev. Lett. 103, 
042301 (2009). 

U. A. Wiedemann, Nucl. Phys. B 588, 303 (2000); 

C. A. Salgado and U. A. Wiedemann, Phys. Rev. D 68, 
014008 (2003). 

A. Majumder, E. Wang and X. N. Wang, Phys. Rev. 
Lett. 99, 152301 (2007). 

V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 
438 (1972); Yu. L. Dokshitzer, Sov. Phys. JETP 46, 641 
(1977); G. Aharelh and G. Parisi, Nucl. Phys. B 126, 
298 (1977). 

E. Braaten and M. H. Thoma, Phys. Rev. D 44, 2625 
(1991). S. Wicks, et ai, Nucl. Phys. A 784, 426 (2007); 
G. Y. Qin, et al., Phys. Rev. Lett. 100, 072301 (2008). 
A. Maiumder. arXiv: 0810.4967. [nucbth]. 
A. Majumder and B. Miiller, Phys. Rev. C 77, 054903 

(2008). 

A. Maiumder, larXiv:0901."45T6 ' [nucbth]. 
A. Majumder, to be published. 

X. F. Guo and X. N. Wang, Phys. Rev. Lett. 85, 3591 
(2000); Nucl. Phys. A 696, 788 (2001). 
A. Majumder, R. J. Fries and B. Miiller, Phys. Rev. G 
77, 065209 (2008). 

J. Casalderrey-Solana, E. V. Shuryak and D. Teaney, J. 
Phys. Conf. Ser. 27, 22 (2005); R. B. Neufeld, Phys. Rev. 
G 79, 054909 (2009). 



